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We analyze the relationship between the mass of a spherical component and the min- 
imum possible thickness of stable stellar disks. This relationship for real galaxies allows 
the lower limit on the dark halo mass to be estimated (the thinner the stable stellar disk 
is, the more massive the dark halo must be). In our analysis, we use both theoretical 
relations and numerical N-body simulations of the dynamical evolution of thin disks in 
the presence of spherical components with different density profiles and different masses. 
We conclude that the theoretical relationship between the thickness of disk galaxies and 
the mass of their spherical components is a lower envelope for the model data points. 
We recommend using this theoretical relationship to estimate the lower limit for the dark 
halo mass in galaxies. The estimate obtained turns out to be weak. Even for the thinnest 
galaxies, the dark halo mass within four exponential disk scale lengths must be more than 
one stellar disk mass. 
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1 Introduction 



A large body of observational data provides evidence for the existence of dark matter, 
which manifests itself only through its gravitational influence on stellar systems of various 
scales (from individual galaxies to clusters of galaxies and superclusters) , determining to 
a great extent their dynamics and structure. The controversy surrounding the following 
two main questions related to the dark matter is still continuing: What is its nature and 
to what extent does it exceed in mass the luminous matter? The first question is far from 
being resolved. As regards the second question, on the scales of individual galaxies, it is 
formulated as follows: Beginning from which regions (inner or outer) does the dark halo 
mass prevail over the luminous mass? 

There are several observational constraints on the mass and extent of the dark halos 
in galaxies; these constraints do not depend on what the nature of the dark matter is. 
In general, they give an estimate of the lower limit for the dark halo mass within a fixed 
radius. For spiral galaxies, this estimate primarily follows from the analysis of the con- 
tributions from various components of the system to its rotation curve for the so-called 
maximum disk model (a classical example of the separation of the contributions from 
the disk and the dark halo to the rotation curve of a galaxy, NGC 3198; van Albada 
et al. 1985). In this model, the dark halo mass within the optical radius of a galaxy 
generally does not exceed the stellar disk mass. When specific systems are investigated, 
a joint interpretation of the data on the rotation curves and radial velocity dispersion 
profiles, along with considerations regarding the marginal stability of stellar disks (see, 
e.g., Bottema and Gerritsen 1997; Khoperskov et al. 2001; Zasov et al. 2001; Khoper- 
skov and Tyurina 2003)^, often raises significantly the lower limit for the dark halo mass 
within a fixed radius (occasionally, interpretation of the data on the stellar velocity dis- 
persion based on numerical simulations increases the dark halo mass by almost an order 
of magnitude compared to its value given by the maximum disk modeP). 

The upper limit for the dark halo mass is much more difficult to constrain. Where the 
region within the optical radius of a galaxy is involved, stability-related considerations are 
invoked. For example, the presence of a regular spiral pattern related to the propagation 
of density waves in galaxies requires a gravitationally "active" stellar disk, which reduces 
the contribution from the spherical component to the total gravitational field of the spiral 
galaxy and to the total mass (Athanassoula et al. 1987). Another important test for the 
presence of a massive halo is the length of the tidal tails in interacting spiral galaxies. In 
interacting systems, we quite often observe extended tidal features stretching to distances 
as large as 40-100 kpc. Based on numerical simulations of interacting disk galaxies with 
dark halos that correspond to models obtained in cosmological calculations (Navarro et al. 
1999), Dubinski et al. (1999) showed that the calculated parameters could be reconciled 
with the parameters of observed systems only for two types of models: models with an 
extended moderate-mass halo and a rotation curve that is determined almost entirely 
by the stellar disk within the optical radius and models with a compact low-mass halo 
that makes the main contribution to the rotation curve in the inner disk. Reconciling 

-'^In these papers, the resuhs of numerical A^-body simulations with a variable disk mass were used to 
explain the observed stellar velocity dispersion for a number of spiral galaxies. 
^The galaxy NGC 891 may be cited as an example (Khoperskov et al. 2001). 
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the rotation curves of specific interacting systems (e.g., NGC 4038/39 and NGC 7252) 
and the sizes of their tidal features with the model parameters leads one to conclude that 
the stellar disks in such systems dominate in mass within the optical radius, which gives 
a strong upper limit on the dark mass. A similar conclusion was also drawn from the 
simulations of the famous Mice galaxy, NGC 4676 (Sotnikova and Reshetnikov 1998). 

Comparison of the various estimates for the dark mass in galaxies shows that the 
upper and lower limits do not overlap (see, e.g., McGaugh and Block 1998). This gives 
reason to revise a number of tests, in particular, the test related to the thickness of the 
stellar disks in spiral galaxies suggested by Zasov et al. (1991, 2002). These authors 
argue that the thickness of a stable stellar disk at a given halo mass is limited below. If 
the disk thickness is smaller than a certain value, then bending instability will lead to 
increasing the vertical velocity dispersion and to a thickening of the system. The growth 
of bending perturbations is stabilized by a massive dark halo: the more massive the halo 
is, the thinner the stable disk can be. As we showed previously (Sotnikova and Rodionov 
2005), not only a massive dark halo has a stabihzing effect on the growth of bending 
perturbations, leaving the stellar disk fairly thin. A low-mass compact bulge also produces 
such an effect. In this case, the lower limits on the dark mass can be mild and consistent 
with the upper limits. In this paper, we analyze the relationship between the stellar disk 
thickness and the dark halo mass that is derived both from theoretical considerations and 
from numerical simulations of the dynamical evolution of thin disks in the presence of 
spherical components with different density profiles and different masses. The derived 
relationship predicts moderate dark halo masses even for the thinnest galaxies. 



2 The stellar disk thickness and the relative mass of 
the spherical component 

When the global parameters of stellar disks are analyzed, the disks are generally assumed 
to be stable against the growth of perturbations in the plane and in the vertical direction. 
The stability conditions for a given dark halo mass limit below the radial and vertical 
stellar velocity dispersions. Since the vertical velocity dispersion is related to the stellar 
disk thickness, it is also limited below at a fixed dark halo mass. 

The disks are commonly assumed to be at the stability boundary. In this case, the 
relative stellar disk thickness and the relative dark halo mass are related to one another. 
This relationship was (probably first) noticed by Zasov et al. (1991) on the following 
grounds. 

For a gravitating layer with a volume density profile that corresponds to an isothermal 
(along the z axis) model (Spitzer 1942), 

p(^) = — sech2(^/^o) , (1) 

the following relation, which is a vertical equilibrium condition, is known to be valid: 

al = ttG^oS . (2) 
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Here, S is the star surface density of the layer, zq is the scale length of the density 
variations in the z direction, and is the vertical stellar velocity dispersion of the layer. 
If condition ^ is applied to a stellar disk, then the R dependences of S, zq, and 
should be taken into account. For real stellar disks embedded in a dark halo, the vertical 
equilibrium condition must be defined slightly differently (Bahcall 1984). However, as was 
shown, for example, by Rodionov and Sotnikova (2006), the approximation of isothermal 
layers holds good for the entire disk except the innermost regions even in the presence of 
a massive spherical component. 

If the disk is stable against perturbations in the plane, then the root-mean-square 
(rms) radial velocity must be 

^ 3.36 GS(i?) 



v'^ V dv 

where (5t~1 is the Toomre parameter, k = -i / 2 ( — ^ H — -—r— ) is the epicyclic frequency. 



RdR, 

and Vc is the linear circular velocity. The coefficient 3.36 in Eq.Q corresponds to an 
infinitely thin model disk (Toomre 1964). 

An expression for the minimum possible relative thickness of a stable disk can be 
derived from equilibrium condition Q and stability condition Q. To be able to compare 
theoretical estimates with observational data, let us write the expression for the relative 
disk thickness at a certain radius, for example, R = 2h, where h is the exponential radial 
disk scale length, which enters into the commonly used law of surface density variations 
with radius: 

S(i?) = ^exp(-i?A), 

where is the total disk mass. The radius 2h was chosen on the grounds that the 
observed vertical disk structure in this region is generally no longer distorted by the 
central bulge, the bar, and the so-called X-shaped structures. 

At R^2h, the rotation curves of luminous galaxies generally reach a plateau. In this 
case, the epicyclic frequency (see stability condition 0) is related to the linear circular 
velocity Vc by k{R) = \/2vc/R- In the plateau region, Vc is roughly constant and, hence, 
Vc{2h) ~ fc(4/i). We can then use Vc to estimate the total mass of a galaxy (including 
the mass of its spherical component Msph) within a fixed radius, for example, R = 4h, 
Mtot(4/i) ~ Ahv'^/G (see, e.g., Zasov et al. 2002). In addition, the marginal stability 
condition in the disk plane allows Qt to be chosen almost uniquely. Based on numerical 
A^-body simulations, Khoperskov et al. (2003) showed that at i? ^ (1 ^ 2)/;, for 
marginally stable 3D stellar disks is equal to a nearly constant value of 1.2 — 1.5, which 
does not depend on the relative mass of the spherical component. Based on the results 
of Khoperskov et al. (2003), we took Q"[{2h) ^ 1.4. Using Eqs.Q and under our 
additional assumptions, we then obtain 

z,{2h) ^ al Md 

h ^ ■ al MU^h) ■ ^ ' 

If we fix the ratio az/a^ at the level given by the local linear criterion for bending 
instability, i.e., at approximately 0.29 — 0.37 (Toomre 1966; Kulsrud et al. 1971; Poly- 
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achenko and Shukhman 1977; Araki 1985), then, according to Eq.Q, the relative stellar 
disk thickness measured at the given distance R = 2h is 



h ^ 1+ fi{4h) ' ^ ' 

where fi{Ah) = Msph(4/i)/Md(4/i), Mtot(4/i) = Md(4/i) + Msph(4/i), Msph(4/i) - mass of 
spheroidal component. In (0), we took into account the fact that 90% of the disk mass is 
contained within Ah for an exponential surface density profile. The numerical coefficient 
in Eq.(jSI) was obtained for az/crn = 0.37 (this value corresponds to the linear criterion 
for bending instability in Polyachenko and Shukhman (1977)). The larger the relative 
mass of the spherical component within a fixed radius, for example. Ah, the smaller 
the relative disk thickness ZQ{2h)/h at this radius. The latter conclusion is one of the 
main conclusions reached by Zasov et al. (1991). In the conclusion under discussion the 
coefficient in relation (0) will be of interest in the subsequent analysis. 

Equation was derived for the specific law of density variations in the vertical 
direction (0). In certain cases, the vertical volume density profile of the disk matter is 
fitted by other laws (van der Kruit 1988), for example, 

p{z) = sech.{z/ze) (6) 

TTZe 

or an exponential law, 

J] 

= ^ exp(-|z|//i^), (7) 

where Ze and h^ are the corresponding vertical scale heights that can depend on R. It can 
be shown that an equilibrium condition similar Q for profiles (0) or ((Tj) may be written 

as 

< >= CttGAS , (8) 

where < o"^ > where is the vertical velocity dispersion averaged along the z axis; A 
is the half-thickness of a homogeneous layer with a density equal to the density in the 
disk at z = 0; and C = 1, C = 28C(3)/7r^ ^ 1.0854 {( is the Riemann function), and 
C = 1.5, respectively^, for density profiles ((l}, (jH)), and ((Zj). Below, we argue that the 
vertical profiles of model disks at late evolutionary stages are more likely described by 
laws (121) or (jni) than by ((Zj). In this case, C ~ 1.0 — 1.1 in equilibrium condition (jHl) and the 
numerical coefficient in (0) is virtually constant if A, the half-thickness of a homogeneous 
layer^, is used in place of Zq. 

Note that relation is valid for a marginally stable stellar disk. If a galaxy has a 
margin of stability against perturbations in the plane or against bending perturbations, 
then the greater sign must be substituted for the sign of approximate equality in As 

^To obtain these coefBcients, we must substitute the corresponding laws of density variations in the 
Jeans equation that describes the vertical equilibrium and, having calculated the law of velocity dispersion 
variations along z, average this quantity. 

^For density profile (O, the vertical scale height zq is exactly equal to the half-thickness of a homoge- 
neous layer A. 
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a result, for stable stellar disks, Eq.® yields 

The lower limit for the mass of the spherical component in a galaxy, in particular, the 
dark halo mass, can be estimated from the observed relative thickness of the galaxy using 
relation Q. Note that the theoretical estimate is weak. Even for extremely thin galaxies 
with zo/h = 1/15, we find that /i(4/i)^1.7. For example, for the thinnest galaxy in the 
sample that was used by Zasov et al. (2002), zo/h ^ 1/11 (see Fig.4b in this paper). 
For such a galaxy, according to ©, /i(4/i)^0.98. In the sample of edge-on galaxies whose 
parameters were discussed by Kregel et al. (2002), z^/h ^ 1/9.7 for the thinnest galaxy 
(435-G25) (see Table 1 in this paper)^. For such a galaxy, Eq.Q yields /i(4/i)^0.75. 



3 Numerical simulations 

The relationship between the relative stellar disk thickness and the relative dark halo 
mass derived from five N-body simulations was also first discussed by Zasov et al. (1991). 
A similar relationship was constructed by Mikhailova et al. (2001) using several tens of 
simulations. In general, the two relationships coincide. 

It should be noted that theoretical relation (0) written above was derived from the 
considerations that were used to substantiate the results of A^-body simulations (Zasov 
et al. 1991, 2002). However, even a superficial comparative analysis shows that the 
curve corresponding to theoretical relation (0) lies well below the curve constructed from 
the experimental data points (see Fig. 2 in Mikhailova et al. (2001)). According to the 
analytical studies of various authors, the factor az/cR appearing in Eq.(0)) lies within 
the range 0.29 — 0.37. We took the maximum value of this factor, 0.37. If we took a 
smaller velocity dispersion ratio, then the differences between the theoretical curve and 
the experimental data points would be even larger. 

The relationship analyzed by Zasov et al. (1991) and Mikhailova et al. (2001) can be 
roughly approximated by the formula (0.3 — 0.35) /(I + /i(4/i)). Thus, the thickness of the 
numerical stellar disk models is approximately twice the thickness calculated using (jSJ. 
In other words, the relationship derived from A^-body simulations predicts a higher dark 
halo mass than estimate In our view, this discrepancy makes it difficult to interpret 
both the results of numerical simulations and the relationships derived from observational 
data (Zasov et al. 2002) and requires a more refined analysis of the numerical simulations 
that primarily includes the stabilizing role of a compact low-mass spherical component 
(Sotnikova and Rodionov 2005). 

3.1 Numerical Models 

Below, we analyze the parameters of several tens of models that we obtained from nu- 
merical A^-body simulations. We used an algorithm based on the hierarchical tree algo- 

^In this paper, exponential profile |7J was used to fit the vertical stellar disk density profile. The 
vertical scale height was defined by hz, which, according to the authors themselves, is related to zq for 
density profile ^ by hz = O.Szq. 
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rithm (Barnes and Hut 1986) and some of the codes from the NEMO software package 
(http:astro.udm.edu/nemo; Teuben 1995). 

The stellar disk was represented by a system of N gravitating bodies {N was varied 
within the range 300 000 to 500 000) and the spherically symmetric component was de- 
scribed in terms of an external static potential that was a superposition of two potentials, 
the bulge and the dark halo. 

The initial particle velocities (the rotational and random velocities) were specified on 
the basis of equilibrium Jeans equations using a standard technique (see, e.g., Hernquist 
1993). A description of the method for constructing the initial model can be found 
in our previous papers (Sotnikova and Rodionov 2003, 2005). Also given there are the 
details concerning the bulge model (a Plummer sphere) and the halo model (a logarithmic 
potential). 

The ratio of the mass of the spherical component to the disk mass within Ah was varied 
between 0.25 and 4. The mass of the spherical component was distributed differently 
between the bulge and the halo. The scale length of the bulge was chosen to be = h/7 
almost for all models. We also considered models with concentrated (ah = h/1.75) and 
"loose" (tth = h/0.35) halos, Oh is the scale length of the mass distribution in the halo. 
The potential softening length e for all models corresponded to our recommendations 
given previously (Rodionov and Sotnikova 2005), e = h/175. 

For all models, we chose nearly equilibrium (in the radial and vertical directions), but 
unstable (against the growth of bending perturbations) initial conditions. During the 
growth of bending instability, the system was dynamically heated, the velocity dispersion 
increased in the z direction, and the disk thickness increased to a value close to the lower 
boundary of a steady state. The main parameter in which the stellar disk models being 
analyzed differ is the fraction of the total mass of the galaxy contained in its spherical 
component (the dark halo and the bulge). 



3.2 Determining the stellar disk thickness in numerical simula- 
tions 

Before discussing in detail the relationship between the relative stellar disk thickness and 
the relative dark halo mass in galaxies, we should gain an understanding of how the 
thickness of numerically simulated disks is determined. For model disks, the thickness at 
a given distance R from the disk center is usually estimated as the rms deviation of the 
z coordinates of stars from the symmetry plane of the disk (see, e.g., Khoperskov and 
Tyurina 2003): 

^rms = ^J{{Z- < Z >f). (10) 

The disk in the z direction is commonly assumed to have volume density profile (0) 
that corresponds to an isothermal model. In this case, Zo and ^rms are related by 

_ 2^3 

Zq — ^rms ~ 1-10 ^rms • (H) 

71 

For density profile (^, zq is the half-thickness of a homogeneous layer. A, with the 
volume density at z = 0. The parameter zo is estimated from observational data for real 
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galaxies by analyzing how the surface brightness of the edge-on stellar disk varies along the 
z axis and by fitting the brightness distribution to cx seQb?{z/ zq). Based on (fTTj) . one can 
assume that the parameter Zy.^s determined from ()10j) for numerical models quantitatively 
differs only slightly from the disk half-thickness zq appearing in law and characterizes 
the thickness of real disk galaxies. However, we will now show that in practice Zrms is 
a poor characteristic of the model disk thickness and that using it can lead to incorrect 
relations between the stellar disk thickness and other global parameters of spiral galaxies. 

A certain number of particles that fly far away from the disk plane appear in a model 
galaxy during its evolution. Since their number is small, the z particle distribution is 
generally well fitted by isothermal density profile (0). However, this small number of dis- 
tant particles contributes significantly to z^^^^ and leads to an overestimation of the stellar 
disk thickness calculated in this way. The extent to which the thickness is overestimated 
depends on the relative mass of the spherical component and on the distance from the 
center at which we measure the thickness. 

Let us introduce another parameter that can characterize the thickness of model stellar 
disks at a given distance R: Zi/2 — the median of the absolute value of z. Twice 2:1/2 is 
the thickness of the layer within which one half of the particles is contained. For density 
profile (PJ, based on the introduced definitions, we have 



A = zo, 

A _ 1 ^ 

22:1/2 In 3 

Z^Yns ^ 



0.91, (12) 
0.83. 



22:1/2 2^3 In 3 

It is easy to show that the relations between zi/2. A, and 2:1.1ns for density profile 
and ((Tj) can be written as follows: for ©, 



A — — ^ 

A TT 



4 In f tan ^ 



0.89 



0.89, (13) 



221/2 



and for (|7j). 



^ = =h 

A 1 

0.72, (14) 



2^1/2 2 In 2 

2^rms 1 



22:1/2 V21n2 



1.02. 



These relations allow us to check how adequately the introduced parameters describe 
the model disk thickness and to reach certain conclusions regarding the law by which the 



8 



1.5 - 



CM 

I — I 

1 



0-5 



2,5 



Nl 



1.5 



0.5 



(a) 



0.91 

Model 60 
Model 76 
Model 75 
Model 49 
Model 53 



7+: 



(b) 



! h 

0.83 

Model 50 

Model 76 - 

Model 75 

Model 49 

Model 53 



I ' 1 1 1 1 ' ^ 

2 4 6 8 10 12 14 

R 



Figure 1: (a) A/(2zi/2) vs. R for several models. The solid line indicates the theoretical 
value of 0.91 for this ratio for density profile (P). All models have the same mass of the 
spherical component and differ only in mass distribution of this component between the 
bulge and the halo (the models are listed in order of increasing bulge mass fraction). For 
all models, Msph(4/i)/Md(4/i) = 0.5 and h = 3.5 kpc. (b) 2;rms/ (22:1/2) vs. R for the same 
models. The solid line indicates the theoretical value of 0.83 for this ratio for density 
profile (HJ. 
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Figure 2: Same as Fig. ^ for two other models - with relatively massive (model 44) and 
very massive (model 12) spherical components. These models are bulgeless. For models 
44 and 12, /i(4/i) = 1.5 and /i(4/i) = 3.43, respectively; h = 3.5 kpc. 
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vertical density profile can be fitted. Figures ^ and El show the azimuthally averaged radial 
profiles of the ratios /S./2zij2 and -Srms/2-2i/2 for several dynamical stellar disk models. 

We see from Figs.^ and|2K that the ratio A/2zi/2 is almost constant along the radius 
for various models. Note that the particle parameters from different parts of the disk 
are used to determine 2:1/2, A, and z-rms- When 2:1/2 is estimated, it does not matter how 
extended the tail of distant particles is, but the number of particles that contribute to Zi/2 
is large and accounts for exactly one half of all particles. The volume density produced 
by only those particles that are located near the main plane of the galaxy and the total 
surface density are used to calculate A; the positions of all the remaining particles do not 
affect in any way the determination of A. As regards 2:rms, all particles, including those 
that randomly flew far away from the disk plane, contribute to this quantity. It follows 
from the constancy of the ratio /S./2zij2 that these two quantities are adequate estimates of 
the model disk thickness; otherwise, it would be difficult to explain the coherent behavior 
of these quantities. At the same time, the equality of this ratio to 0.9 is an argument for 
the fact that the vertical density proflle of stellar disks is similar to law (P) or (jHl). 

As follows from Figs. HJd andEb, ^rms/2zi/2 varies signiflcantly along the radius R for 
all of the analyzed models and, in addition, this ratio is almost always systematically 
higher than that given by theoretical relations (fT^ . p3|) . and even (fT^ . The thickness 
overestimation strongly depends on the degree of matter concentration in the spherical 
component and on its relative mass. It is different in different parts of the disk. The 
differences are largest for the models with bulges. For them, the disk thickness in a region 
of the order of two exponential disk scale lengths can be overestimated by 100% or more 
(models 49 and 53, Fig.^) if ^rms is used. As regards regions of more than two radial disk 
scale lengths, the thickness is overestimated here by 50% or more. For the models with a 
relative mass of the spherical component iJ,{Ah) < 1.5 without a bulge (models 50 and 44, 
Figs. ^ andEb), the thickness of the central regions is more or less adequately described 
by the parameter z^ms- In this case, the density proflle corresponds to isothermal law ((H), 
for which 2;rms/22i/2 ~ 0.8. However, at distances larger than one or one and a half 
radial exponential disk scale lengths, the thickness proves to be grossly overestimated, 
approximately by 50%. Only for the models with /i(4/i)^3.0 (model 12, Fig. ^) does 
Zt-^s serve as a good estimate of the thickness. In addition, the Zj-ms fluctuations along R 
for the model disks are too large even when a large number of particles is used, in contrast 
to 2:1/2 (Sotnikova and Rodionov 2005). The ratio z^^^jlzxj^ also exhibits the same large 
fluctuations. 

The parameter z^ is commonly used to analyze observational data by assuming the 
vertical density proflle to be isothermal. For our numerical models, we can determine zq 
via both z^.^^ (zq ~ 1-10 -Zrms) and ^1/2 (2:0 ~ 1.82 ^1/2). As we see from the plots of the 
radial disk thickness proflle (Figs.Eti' andEfc), 2:9 calculated considerably larger 

than 2:0 calculated via 2:1/2. Figures and Eli show the vertical density proflles at a given 
radius {R = 3.5 kpc, which corresponds to the exponential radial disk scale length) and 
the flt to the density proflle by isothermal law (0) with the parameter zq calculated as 
zq ^ 1.10 2rms and as zq ~ 1.822:1/2. We see that if zq is determined via 2:1/2, then law ((H) 
flts excellently the vertical density proflle in the model. When zq is determined via Zj-ras, 
the discrepancies between the proflles are enormous. 

The above arguments lead us to the following conclusion: the parameter 2:1/2 should be 
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Figure 3: Comparison of the disk thicknesses calculated using 2;rms and 2:1/2 for two models 
at the time t = 5000 Myr. Models 49 and 53 have M,ph(4/i)/Md(4/i) = 0.5, where h = 3.5 
kpc, but differ in bulge mass fraction (it is larger for model 53). (a), (c) Zq vs. R; the solid 
and dashed lines represent zq estimated as 1.82^1/2 and 1.10 2;i.ms, respectively, (b), (d) 
Vertical particle density profiles at the radius R = 3.5 kpc; the thick solid line represents 
the profile obtained from the model; the thin solid line represents profile in which zq 
was estimated via zi/2] the dashed line represents profile in which zq was estimated 
via Zrms. The disk surface density at a given distance R in Eq. (jT)) was calculated from 
the model. 
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taken as the stellar disk thickness to analyze the relationships between the disk thickness 
and other global parameters of galaxies in numerical models. Using Zj-ms can distort 
significantly the relationship between the relative thickness of a marginally stable stellar 
disk and the relative mass of the spherical component and can raise the lower limit for the 
dark halo mass. Below, we use 1.82 2:1/2, since it corresponds to zo for density profile (jTJ. 

3.3 Analysis of numerical simulation results 

In deriving Eq. for the ratio of the vertical and radial velocity dispersions, we used the 
value given by the local linear criterion (Toomre 1966; Kulsrud et al. 1971; Polyachenko 
and Shukhman 1977; Araki 1985). However, the saturation level of the bending instability 
in numerical simulations can be considerably higher than this value (see, e.g., Sellwood and 
Merritt 1994; Merritt and Sellwood 1994; Sotnikova and Rodionov 2003). This probably 
explains why the relationship in Zasov et al.(1991) and Mikhailova et al.(2001) lies above 
the "theoretical" one. 

The ratio cXz/cTfj at which the bending instability is saturated depends primarily on 
whether the system admits the growth of global modes. If the growth of global modes is 
suppressed, then the ratio (Tz/ctr at a given radius is determined by the local conditions 
and is equal to its value that follows from the linear criterion. If, however, global modes 
are possible in the system, then the level of vertical disk heating can be considerably 
higher. 

Sweeping throughout the disk, the global modes cause the ratio <Jz/<Jr to increase to 
0.6 and, in central disk regions, to 0.7 — 0.9 (Sotnikova and Rodionov 2003). In Fig. |3J 
the ratio OzjoR is plotted against the maximum amplitude of bending modes in the 
time interval from to 3000 Myr for the distance R = 2h. The plot was constructed 
using numerical simulations of disk galaxies. We see a clear correlation between the 
maximum amplitude of the growing and propagating (through the disk) perturbation and 
the instability saturation level and, accordingly, the final az{2h)/aii{2h). When global 
modes are possible in the system, the maximum perturbation amplitude turns out to be 
larger. 

For inhomogeneous disks, the possibility of the growth and propagation of global modes 
is determined by the conditions in the central disk regions (Sellwood 1996; Sotnikova and 
Rodionov 2005). These conditions impose constraints on the vertical oscillation frequen- 
cies related to the spherical component and the disk and for the model with /i(4/i)^2.0 
depend not only on fi, but also on the degree of matter concentration in the spherical 
component toward the center (Sotnikova and Rodionov 2005). A low-mass compact bulge 
suppresses the growth of global bending modes as effectively as a massive dark halo with 
a large scale length of the matter distribution. 

In Fig. EK, o'z{'^h)/aji{2h) is plotted against the relative mass of the spherical compo- 
nent iJ,{Ah) for a stellar disk that came to equilibrium after the saturation of the bending 
instability {t = 5000 Myr). It should be noted that a difficulty with the determination of 
the exponential stellar disk scale length h arises when constructing this plot. During the 
evolution, an appreciable mass redistribution takes place in many models. The density 
in the central region can increase greatly. The most significant changes are observed in 
the models with the formation of a bar. As a result, the density profile in these central 
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Figure 4: az{2h)/aji{2h) determined at the time t = 3000 Myr vs. maximum disk warp 
amplitude A^^ax for the period from to 3000 Myr (the disk region within 10 kpc was 
considered) . 

regions was fitted by an exponential law with a considerably smaller scale length than 
that of the density profile on the periphery. For our analysis, we took the exponential 
scale length estimated from the periphery of a model disk (the region from 5 to 10 kpc), 
Fig. EK- For comparison. Fig. shows the same plot, but the value of h was taken at 
the initial time {h = 3.5 kpc for all models at the initial time). The total mass within 
4:h also varies with h. This explains the differences in the ranges of /i(4/i) for the same 
set of models in Figs. and Eb- The squares in both figures indicate the data for the 
models with a bulge whose mass is several times lower than the disk mass. Irrespective of 
/i(4/i), the minimum value of az{2h) /(TFi{2h) obtained from the models satisfies the linear 
criterion for bending instability, (0.29 — 0.37). This minimum value corresponds to the 
models with a compact bulge or a massive halo. 

The stabilizing role of a compact bulge lies in the fact that it inhibits the growth 
of large-scale high amplitude perturbations. In the models with a low mass spherical 
component and without a compact bulge, the processes that may be arbitrarily called a 
"violent" stage of the bending instability proceed during the evolution (for the appearance 
of a bell mode in hot models with large Qt and the rapid growth of bending instability of 
the bar in cold disks, see Sotnikova and Rodionov (2003)). The physics of these processes 
may not be directly related to the physics of the bending instability. A vertical resonance 
can appear in such models (Combes and Sanders 1981; Combes et al.l990; Pfenniger and 
Friedli 1991; Patsis et al.2002). Nevertheless, these processes affect <Jz/<Jr — it becomes 
considerably higher than that given by the linear criterion for bending instability. 

In Fig.inti', ZQ{2h)/h is plotted against yu(4/i) at the time t = 5000 Myr for h calculated 
from the region from 5 to 10 kpc. Figure shows the same plot for an initial value 
of h = 3.5 kpc. As expected, at fixed relative mass of the spherical component fi{4:h), 
the models with bulges can be considerably thinner than the models without bulges 
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Figure 5: az{2h)/aji{2h) vs. /x(4/i) for our models at the time t = 5000 Myr. The 
crosses and squares indicate the models without and with a bulge, respectively. The 
dashed line indicates OzjoR — 0.37 that corresponds to the saturation level of the bending 
instability derived from the linear criterion of Polyachenko and Shukhman (1977). (a) The 
exponential disk scale length h was calculated from the model by the least-squares method; 
the disk region of > 5 kpc and < 10 kpc was used to calculate h. (b) The initial 
exponential disk scale length h — 3.5 kpc was taken as h. 
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(see Figs. and IHb)- This effect is most pronounced for the models with a spherical 
component of low relative mass (/i(4/i)^l). 

The following result proved to be more unexpected. The analytical relation between 
the minimum possible zo{2h)/h for a stable disk and /i(4/i) (see 0) is a lower envelope 
for the model data (see Figs. IHK andlHb)- In other words, at fixed /i(4/i), ZQ{2h)/h for the 
thinnest model is approximately equal to the value obtained from Eq.©. 

The data points lying on the lower envelope correspond to the models with compact 
bulges or with massive halos. For such models, the bend was formed in "quiet" regime 
or was suppressed altogether. In this case, the final value of (Jz/crR was also close to the 
saturation level of the bending instability derived from the linear criterion. If, however, 
the bending instability in the model grew in "violent" regime, then the final value of 
<^z/<^R could be considerably higher than that given by the linear criterion and, hence, 
the final thickness of the model galaxy could be much larger than that following from the 
linear criterion. 

Thus, if a compact spherical component that can be represented by both a bulge and 
a dark halo with a sharp central density peak (a cuspy halo) is present in a real spiral 
galaxy in its central region, then the bending instability in such a galaxy will grow in quiet 
regime. In this case, if the galaxy is marginally stable against bending perturbations and 
perturbations in the pane, then its parameters will obey Eq. (0). 

Consequently, we can use only relation (jH]) to estimate the lower limit for the dark halo 
mass from the relative thickness of galaxies. We emphasize once again that this relation 
gives a weak constraint on the dark halo mass, and only on the condition that the galaxy 
is in equilibrium, even for thin galaxies. 

4 Discussion and conclusions 

It is well known that for a stellar disk to be stable against perturbations in the plane, 
the value of (Tr must be higher than a certain critical value (Toomre 1964). On the other 
hand, the stability against the growth of bending perturbations requires that the ratio 
Cz/cr be also larger than a certain critical value. Consequently, the value of is limited 
below for a stable stellar disk. The presence of a spherical component, for example, a dark 
halo, has a stabilizing effect and reduces the minimum radial velocity dispersion required 
for the stability against the growth of perturbations in the plane. It thus follows that 
the disks embedded in a massive dark halo can have lower values of and, as a result, 
can be thinner. These considerations underlie the conclusion about the existence of a 
relationship between the relative thickness of marginally stable stellar disks, z^/h, and 
the relative mass of the spherical component, /i, within a fixed radius (Zasov et al.l991, 
2002). 

In this paper, we estimated the coefficient in this relationship and wrote the corre- 
sponding theoretical relation (jSj). The relationship was confirmed by N-body simulations 
of disk galaxies. This theoretical relation can be used to estimate the lower limit for the 
dark halo mass from the observed thickness of real spiral galaxies 0. However, even for 
the thinnest galaxies, this relation yields a weak constraint on the dark halo mass (the 
dark halo mass within four exponential disk scale lengths must be higher than approx- 
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Figure 6: Relative stellar disk thickness at radius 2h vs. relative mass of the spherical 
component within 4h for our models at the time t = 5000 Myr. The value of zq was 
calculated as 1.82 The crosses and squares indicate the models without and with 
a bulge, respectively. The dashed hne represents theoretical relation (0). (a) The expo- 
nential disk scale length h was determined from the model by the least-squares method; 
the disk region of R > 5 kpc and R < 10 kpc was used to calculate h. (b) The initial 
exponential disk scale length h = 3.5 kpc was taken as h. 
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imately one disk mass). We can adduce several arguments that will relax this already 
weak constraint even further. 

First, we formally obtain an estimate of the lower limit on the total mass of the galactic 
spherical component from Eq. (jHI). The dark halo gives only a partial contribution to this 
mass. Thus, the dark halo mass estimate becomes even weaker. Second, as was already 
noted above, in deriving Eq. (jSl), we assumed that Oz/cyR = 0.37, which corresponds to 
the local linear criterion for bending instability obtained by Polyachenko and Shukhman 
(1977). However, slightly differing values of this ratio, within the range 0.29 to 0.37, were 
obtained in different papers in which the bending instability criterion was derived from a 
local linear analysis. We used the maximum theoretical value for the dispersion ratio. If 
a smaller ratio is taken, then the coefficient in Eq. (jH)) will decrease and, as a result, the 
theoretical estimate of the dark halo mass will also decrease. 

Zasov et al. (1991) and Mikhailova et al. (2001) also derived the relationship being 
discussed from iV-body simulations of disk galaxies (in general, the relationships given in 
these papers coincide). As was already noted above, the model data points in these papers 
(in contrast to the results of our simulations) lie well above relation ^ and, as a result, 
the model relation at fixed stellar disk thickness predicts a considerably higher mass of 
the spherical component than the theoretical relation. What causes these differences? We 
see three causes. 

The main cause of the differences is that the saturation level of the bending instability 
in simulations can be much higher than the level given by the linear criterion. This is 
probably the reason why the relationship being discussed derived in the above papers 
lies above that found from the linear stability criterion for bending perturbations. The 
numerical simulations by Zasov et al. (1991) and Mikhailova et al. (2001) disregarded 
the fact that a compact, but not necessarily massive, spherical component could effec- 
tively suppress the bending instability. As a result, at fixed total mass of the spherical 
component, the galaxies with a more concentrated spherical component toward the center 
(e.g., if there is a compact bulge or a cuspy dark halo in the galaxy) can be considerably 
thinner than the galaxies with a looser spherical component. 

In addition, the overestimation of the model disk thickness is related to the incorrect 
technique of determining the vertical scale height used in the mentioned papers. It follows 
from our analysis that the characteristic disk thickness zq in numerical simulations must 
be estimated via the median of the absolute value of z (-21/2), not via the rms deviation 
of the z coorrdinates of particles from the symmetry plane of the disk (^rms)- Otherwise, 
the disk thickness can be overestimated significantly, by 100% or more. Thus, the mass 
of the spherical component determined from the thickness of a marginally stable stellar 
disk will also be overestimated. Note also that in the mentioned papers, the thickness 
averaged over the entire disk, including its central regions, was taken when constructing 
the experimental relationship. This disk-averaged thickness is definitely larger than that 
calculated on the disk periphery, for example, at the distance R = 2h (for the radial 
profiles of the final model disk thickness, see Mikhailova et al. (2001) and Sotnikova and 
Rodionov (2003, 2005)). 

The ambiguity of the relationship between the thickness of a galaxy and the mass of its 
spherical component with different density profiles in numerical simulations complicates 
the estimation of the dark halo mass from the relative stellar disk thickness Zq / h. However, 
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to estimate the lower limit for the dark halo mass from the relative thickness of real 
galaxies, we can use theoretical relation Q, which coincides with the lower envelope of 
the numerical relationship (Figs. andlHt")- This relation gives a weak constraint on 
the dark halo mass, n{4h)^l, even for thin galaxies, which is consistent with the upper 
limits that follow from considerations concerning the conditions for the existence of spiral 
density waves in galaxies (Athanassoula et al. 1987) and from arguments regarding the 
extent of tidal features in interacting galaxies (Dubinski et al. 1999). 
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